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Abstract 

We present a Next-to-Leading Order calculation of the cross section for the 
leptoproduction of large- E±_ hadrons and we compare our predictions with HI 
data on the forward production of 7r°. We find large higher order corrections 
and an important sensitivity to the renormalization and factorization scales. 
These large corrections are shown to arise in part from BFKL-like diagrams 
at the lowest order. 



1 Introduction 



The electroproduction of \arge-E± hadrons which is observed by the HERA 
experiments HI [TJ |21 E] and ZEUS |3] may provide important tests of QCD. In 
addition to the study of the partonic subprocesses, of the parton distributions in the 
proton and of the fragmentation functions, it also offers the possibility to observe the 
virtual photon structure function. A contribution of the latter is indeed expected 
when the hadron transverse energy squared E\ is much larger than the photon 
virtuality Q 2 = \q 2 \ ; in this case the virtual photon structure function contribution, 

E 2 

proportional to log -pp-, can be important. 

Another interest of this reaction is the study of the production mechanisms of 
forward hadrons. Indeed, the forward region can be associated with BFKL dynam- 
ics and provide tests of the latter j^j. Several papers have studied the production 
of large- E± jets and hadrons [Zj in the forward direction, and have concluded 
that these reactions are relevant for the study of BFKL dynamics. However, results 
from HI P] show that theoretical predictions based on the DGLAP dynamics, im- 
plemented in RAPGAP jHj or based on BFKL dynamics pj are both in agreement 
with data. These theoretical results do not necessarily contradict each other, since 
the same Feynman graphs may contribute to both of them. But it is clear that a 
deeper understanding of the underlying forward dynamics requires quantitative pre- 
dictions, which are not fully realized in the existing literature above. Ref. [B] which 
implements the DGLAP dynamics rests on a Leading Order (LO) approximation 
and suffers from scale dependence which forbids an absolute normalization. Ref. [Oj 
takes into account some Higher Order (HO) corrections to the LO BFKL equation, 
but does not include non-BFKL contributions. 

In this paper we calculate the HO corrections to the Born subprocesses associated 
with the electroproduction of large- E± hadrons namely the QCD Compton process 
7* + q — > g + q and the fusion process 7* + g — > q + q. The Born subprocess cross 
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sections, of order 0(aa s ), and the HO subprocess cross sections, of order (D(aa 2 ), 
are convoluted with parton distributions and fragmentation functions calculated at 
the NLO approximation. The total NLO cross-section (Born + HO contributions) 
is sensitive to the choice of the renormalization and factorization scales, but there 
is a compensation between the variations of the Born and HO contributions in such 
a way that the NLO cross section is more stable than the LO (containing only the 
Born terms) cross sections. 

The work presented here is a fixed order (for the subprocess cross section), NLO 
calculation, the DGLAP dynamics being included by the scale dependent distri- 
butions and fragmentation functions. However, among the HO contributions, one 
of them corresponds to the lowest order BFKL cross section, namely, the reaction 
7* + g — > g + q + q. As this contribution is part of an HO calculation, we have a way 
to establish a link between the normalization of our NLO cross section and that of 
the (lowest order) BFKL cross section. 

As mentioned in the beginning of this section, the HO calculation generates a 
contribution proportional to the virtual photon structure function. At order 0(aa 2 ) 
we obtain the Born expression of this structure function, proportional to log E\jQ 2 
and we shall study how important this contribution to the large- E± forward hadron 
cross section is in the HI kinematical configuration. All order contributions to 
the virtual photon structure function can be resummed using an inhomogeneous 
DGLAP evolution equation. Here we shall briefly discuss this possibility, leaving for 
another publication [TU] a detailed analysis of the resummed virtual photon structure 
function. 

As a final point concerning the nature of this calculation, let us emphasize the fact 
that it describes the production of large- E± hadrons in the virtual-photon proton 
centre-of-mass system (CMS) and that it is also valid in the limit Q 2 = 0, the large 
scale then being provided by E\. It must be compared with experimental results 
which impose a lower bound on the final hadron E± in the 7* — P CMS, as is done 
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by the HI and ZEUS collaboration O El H] • Therefore the present work does not 
consider the target fragmentation mechanism, which requires the introduction of 
fracture functions ^ ^ s a ^ so different from the inclusive calculations 

of ref. f3J E3] m that it is an exclusive NLO partonic generator. This allows us to 
calculate various types of correlations (for example between large- E± hadron and 
jets) and facilitates the implementation of experimental cuts. 

In the next section we shall present an overview of the relevant DIS kinematics 
and of the method used to calculate the HO corrections. Section 3 is devoted to a 
discussion of the virtual photon structure function and, in section 4, we compare our 
theoretical results with HI data [2|. We shall discuss in detail the importance of the 
virtual photon contribution and of the BFKL-like contribution in the HI kinematical 
region. Section 5 studies the production of large- E± hadrons in the central region 
in rapidity. Section 6 is the conclusion. 

2 The NLO calculation 

The kinematics of the reaction e(£) + p(P) — > e(£') + h(P±) + X is fixed by the 
observation, in the laboratory frame, of the outgoing lepton and hadron h momenta 
We define the photon variables (axis Oz along the initial proton momentum) 



q 2 = - q 2 = -{i - ay 

_ gH _ q° ~q z _ P-q _Q 2 1 () 

y ~ £(-) £0_£z p. £ S x Bj ' U 

where we have neglected the proton mass and used the notation S = (£ + P) 2 and 
%Bj = ^~ q - The outgoing hadron is defined by its transverse energy E^2 b and its 
pseudo rapidity r}\ ah . 

The inclusive cross section is written in terms of the leptonic tensor (summed 
and averaged over the spins of the leptons) 
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£ HU = 2 Ufigv + ^ t _ gii v U.g _ m 2^ ; (2) 
and the hadronic tensor T MZ , which describes the photon-proton collision 



d ° " 1 1 rl ^-dPS, (3) 



d(pdQ 2 dy 2vr 2vr 25 J 2 g 4 



where 



dPS = (2tt)V [q + P-fp^ ft 

is the final state hadron phase space element and (p the photon azimuthal angle. (A 
sum over the number of final hadrons is understood in (JSJ)). 

The hadronic tensor can be calculated as a convolution between the partonic ten- 
sor which describes the interaction between the virtual photon and the parton of 
the proton, and the parton distribution in the proton G a (x, M). The fragmentation 
of the final parton which produces a large- E± hadron is described by the fragmen- 
tation function D£(z, Mp). These distributions depend on the factorization scales 
M and M F , 



J T^dPS = £/ ^G a (x, M) J dz D h h {z, M F )tf„ ■ dps, (4) 

a, b 

where dps is the phase space element of the partons produced in the hard photon- 
parton collision. From expressions Q and (J1J), we obtain 



da E^f a ^ r . „ , . r dz 



£ J dxG a (x,M) J -^D h b {z,M B 



dipdQ 2 dydElfdril ab ~ 2n 2tt 

ab 

^dps', (5) 



d^ ab 1 1 i^tf 

2tt (Att) 2 2xS g 4 

where the phase space dps' no longer contains parton 4 which fragments into h(P±). 
Up to this point we have been writing E±± and n]\ ah to emphasize the frame in 
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(a) (b) 



Figure 1: The QCD Compton subprocess (Fig. la) and real (virtual) HO corrections 
to it (Fig. lb (Fig. lc)). 

which the reaction is observed. Of course (JHJ is valid in any frame and from now on 
we shall drop the index "Lab" . 

The tensor product is a series in a s . Taking into account the first and second 
order contributions, we rewrite (jHJ) as 



The cross sections a^ >rn are the subprocess Born cross sections which describe the 
electroproduction of a large- E± parton b and K^ b ° are the associated Higher Order 
corrections. 

Figure la shows the Born term corresponding to the QCD Compton (QCDC) 
process (here partons a and b are quarks). When a is a gluon, we have the so called 
photon-gluon fusion Born term. Examples of a graph contributing to HO corrections 
to the QCDC term are shown in Fig. lb, c. (The numbers 1 to 5 label the initial 
and final partons according to a convention used in the HO calculations described 
below). 

In expression (JHJ) we have explicitly written the dependence of the cross section 
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on the strong coupling constant a s {p) which depends on the renormalization scale \x. 

It is convenient to perform the calculation of a Born and K in the virtual 
photon-proton center of mass system, and from now on we shall work in this frame. 
In fact, the HI collaboration explicitly uses this frame to place cuts on the outgoing 
hadron transverse momentum. We take the positive Oz axis along the proton mo- 
mentum (as per the HI convention) and the leptons are contained in the Oxz-pl&ne. 

It is instructive to give a more explicit form to the tensor product in the 7* — p 
frame using expression (J2J) and defining the transverse polarization vectors = 
(0,1,0,0), e 2 = (0,0,1,0) and the scalar polarization vector = — h^(q z , 0, 0, q°) 
with = (q°, 0, 0, q z ) being the virtual photon momentum 



= Q'itn + t 22 ) + 4 ( Q {1 y2 V) - mlj t n 

+4— 4v^ t 8 i + Q 2 A{1 ~ y) t ss , (7) 
y v y 2 

where y = foZ% — ^3 is identical to the variable defined in (JTJ) in the Lab frame. 

In the limit Q 2 — and after azimuthal averaging over ^4 we recover the unin- 
tegrated Weizsacker- Williams expression 

1 = ( l + (l-y) 2 _ 2y m 2 e 

2 { yQ 2 

with cr ± = i (t n +t 22 ). 

Actually the limit (jHJ) is correct only if lim t ss = 0(Q 2 ). This is not true if an 
initial collinearity is present in the partonic tensor (light partons are massless) which 
leads to the behavior lim t ss = 0(1). This point will be discussed in a forthcoming 
publication [TU] . 

After these kinematical preliminaries, let us describe the calculation of the HO 
corrections which uses the phase-space slicing method elaborated in ref. ^H]- We 
outline the strategy only briefly; for more details we refer to JH] • 
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For a generic reaction 1 + 2^3 + 4 + 5 (fig. 1), at least two particles of the final 
state, say 3 and 4, have a high E± and are well separated in phase space, while the 
last one, say 5, can be soft, and/or collinear to any of the four others. Of course on 
the photon side there is no collinear divergence as long as Q 2 is different from zero. 
This part of the calculation is related to the virtual photon structure function to be 
discussed below. In order to extract the other singularities, the phase space is cut 
into two regions : 

(1) Part I where the norm E±$ of the transverse momentum of particle 5 is 
required to be less than some arbitrary value E± m taken to be small compared to the 
other transverse momenta. This cylinder contains the infrared and the initial state 
collinear singularities. It also contains a small fraction of the final state collinear 
singularities. 

(2) Parts Ila(b) where the transverse momentum vector of particle 5 is required 
to have a norm larger than E± m , and to belong to a cone C^{C^) about the direction 
of particle 3(4), defined by (775— ?7i) 2 + (05 — (pi) 2 < Rfh = 3, 4), with R t t some small 
arbitrary number. C${C±) contains the final state collinear singularities appearing 
when 5 is collinear to 3(4). 

(3) Part lie where E±§ is required to have a norm larger than E± m , and to belong 
to neither of the two cones C3, C4. This slice yields no divergence, and can thus be 
treated directly in 4 dimensions. For this regular part of the calculations we use the 
cross sections from [T7j . 

The contributions from regions I and Ha, b are calculated analytically in d = 
4 — 2e dimensions and then combined with the corresponding virtual corrections 
(borrowed from DISENT [Hj) such that the infrared singularities cancel each other, 
leaving only the initial (final) state collinear singularities, which are factorized and 
absorbed into the parton distribution (fragmentation) functions. The MS factor- 
ization and renormalization schemes are used in this calculation. 

After the cancellation, the finite remainders of the soft and collinear contributions 
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in parts I and Ha, b, c separately depend on large logarithms In E± m , In 2 E± m and 
In R t h- When combining the different parts, the cancellations of the E± m and 
Rth dependent terms occur. Actually, in part I, the finite terms are approximated 
by collecting all the terms depending logarithmically on E^ m and neglecting the 
terms proportional to powers of E± m . Similarly in parts Ha and lib we keep only 
the logarithmic terms In R t h- Therefore the parameter E± m must be chosen small 
enough with respect to E ±4: so that the neglected terms can be safely dropped. On 
the other hand, it cannot be chosen too small for then numerical instabilities may 
occur. Similar remarks are also valid for the Rth cut-off. 

This approach allows us to build a partonic event generation which is very flexi- 
ble ; various sorts of observables and experimental cuts being easily handled. More 
references to this method, which has been used to calculate NLO correction to sev- 
eral photoproduction and hadroproduction reactions can be found in ref. [THj . 

3 The resolved contribution 

As mentioned in the Introduction, the calculation of the HO corrections leads to 

E 2 

a contribution proportional to log -M- (when Ej_ A 3> Q 2 ) : the so-called resolved pho- 
ton contribution. Indeed, let us consider the contribution associated with Fig. lb in 
which we interchange the label 5 and 4. The integration over E±s (the unobserved 
final quark momentum) produces, among other contributions, a logarithm contri- 

E 2 

bution log -M- associated with a configuration in which the final quark is collinear 
to the virtual photon. More explicitly, we obtain the following expression (for a 
transversely polarized photon) 



0± 









/ dz 


'z 2 + (1 - zf 


ft 









a 
2^ 



+z(l - z)q 1 



^4 



dk\ 



E 2 . 



dk\ 



(k\ - q 2 z{l - z)) 



q 2 z(l — z) 
2 - \ a(0)dps, 



(9) 



S 



in which we have defined k = q — p§ and z = -j=y. The cross section of the 2 — > 2 
subprocess is calculated with k on-shell (k 2 = 0) and dps is the final partonic phase 
space of parton 3 and 4. This expression, as well as a similar one for the scalar 
cross section, will be derived in ref.^U]. A discussion of the second term of Q (the 
non-logarithmic piece) is also postponed to this paper. Here we are interested in 
the term proportional to log E\JQ 2 and in a discussion of the upper limit of the 
integral in (JHJ). 

From Q we obtain the term 



(with the definition P qi {z) = ^{z 2 + (1 — z) 2 ) which defines the quark distribution 
in the virtual photon 



In expression (jllj) . the quark distribution is calculated with no QCD correction. 
However when E\± 3> Q 2 it becomes important to calculate these corrections and to 
replace (jllj) by the LO or NLO expressions of the quark distribution [TTH lT9~j l2T)| l2*T] . 

In this paper we are only interested in the study of the importance of the resolved 
contribution obtained in the kinematical configuration of the HI experiment and we 
content ourselves with the lowest order expression (jllj) . However this expression, 
obtained in the limit Ej_ 4 ^> Q 2 , must be elaborated in order to cover other kinemat- 
ical configurations as well. When Q 2 > E]_ A , we cannot neglect the k\ dependence 
of a{k 2 ± )dps which suppresses the logarithmic integration in expression (jnj). The 
resolved contribution coming from the collinear configuration can be approximated 
by the form 




(10) 



q^z.E^Q 2 ) = P qi (z) log 



Q 2 



(ii) 




(12) 
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which has the correct limit for E\ A ^> Q 2 and Ej_ 4 < Q 2 . In this work we define 
the resolved contribution by the expression (|12j). and in Section 4, we shall calculate 
its numerical importance. 

When the lowest order expression ()12j) is resummed, it must first be removed 
from the HO corrections. We shall call HO s the remaining corrections, and we shall 
say that the subtraction has been performed at the scale M 2 = Q 2 + E 2 LA . 

4 Results 

In this section we present the results obtained with the NLO code described in 
sections 2 and 3. We shall compare our predictions with a selected set of HI data 
|2j and we shall concentrate on a detailed discussion of the various contributions to 
the cross section, namely the HO corrections, the virtual photon structure function 
contribution and the BFKL-like contribution. Here we do not intend to perform 
a complete phenomenological study of the HI data 0] that we shall present in 
a forthcoming publication. We use the MRST99 (higher gluon) distribution 
functions corresponding to = 300 MeV and the KKP fragmentation functions 
of quarks and gluons in ir° j2H]. The renormalisation and factorization scales are 



First we study the cross section da/dxBj measured by HI [2] in the range 
4.5 GeV 2 < Q 2 < 15 GeV 2 with a lower bound on the transverse energy, in the 
7* — p frame, of the forward 7r° given by E±4 > 2.5 GeV. The HERA proton and 
electron beams have laboratory energies 820 GeV and 27.5 GeV respectively, and 
the inelasticity defined in (0) is restricted to the range .1 < y < .6. The for- 
ward domain in which the meson is observed is given by 5° < 6% ab < 25° and 
x,,- = E^ ab j Ep,^ ton > .01. Then we shall consider other kinematical ranges for Q 2 . 

The HI data are compared to our predictions in Fig. 2. The HO contributions 
from which the resolved contribution has been subtracted (see the discussion in 



taken equal to JQ 2 + E\ 
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Figure 2: Comparison with HI data in the range 4.5 GeV < Q < 15 GeV . 
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Figure 3: Comparisons with HI data for two Q 2 -ranges : 2 < Q 2 < 4.5 GeV 2 (left) 
and 15 < Q 2 < 70 GeV 2 (right). 

section 3) are indicated by HO s . In Fig. 2 the importance of the HO corrections, 
especially at small values of Xsj is dramatically visible. As discussed below, it is 
associated with the opening of new channels. 

The resolved component is non- negligible in this range of E\ A and Q 2 . Of course 
its amplitude depends on the factorization scale, which we discussed in Section 3. 
We recall that the factorization scale used here is M 2 = Q 2 + Ej_ 4 . It is interesting 
to look at < Ej_ 4 > and compare it with the value of < Q 2 >. According to 
our calculation < Ej_ 4 >~ 15.3 GeV 2 , which does not depend on the value of 
x B j', < Q 2 >= 6.15 GeV 2 and 10.36 GeV 2 for the ranges 1.10" 4 < x Bj < 2.10" 4 
and 5.5. 10~ 4 < XBj < H-0 10~ 4 respectively. These values lead to a virtual photon 

— ) — 1.25, and 0.91. The dependence 
of < Q 2 > on Xbj also explains the relative decrease of the resolved component at 
large x B j- For two other Q 2 ranges we compare HI data with theory in Fig. 3. 
We clearly notice the decrease of the resolved component when Q 2 increases and 
becomes larger than Ej_ 4 . 
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Figure 4: Examples of HO diagrams which contribute to the "BFKL Born" term 
when the detected hadron is a fragment of the gluon or of the quark q. 

Let us end this general discussion of our results by noting that theory under- 
estimates data by a small amount. Nevertheless we must keep in mind that two 
points are still missing for a more complete comparison. First we have not yet 
studied the scale dependence of our results which have been obtained for the choice 
fj, = M = Mp = M 7 = (Q 2 + Ej_ 4 )^, and second, we have not considered HO 
corrections to the resolved contribution which are known to be large. Indeed we 
can estimate these corrections by using the Weizsacker- Williams approximation (jSJ) 
which is implemented in the photoproduction EPHOX code^S]- From this code we 
obtain a ratio HO /Born ~ 1 in the kinematical domain corresponding to Fig. 2. 

Now we turn to a detailed study of the HO contributions. As expected the two 
contributions corresponding to the Feynman graphs shown in Fig. 4 are the largest 
in the forward direction, because of the exchange of a gluon in the t-channel. They 
also correspond to new subprocesses that are not present at the Born level, as soon 
as the observed partons are the final gluon or quark q ; therefore they do not possess 
singular collinear configurations of partons which contribute to the dressing of the 
distribution and fragmentation functions already present in the Born terms. These 
graphs, with a trigger on the gluon or quark q, also correspond to the Born terms of 
the BFKL ladder in which extra gluons are emitted by the t-channel gluon. This is 
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Figure 5: Comparison of BFKL Born contributions with the total NLO cross section 
and HI data for 4.5 GeV 2 < Q 2 < 15.0 GeV 2 . 



precisely the contribution to the forward cross section [JJ |S] that HERA experiments 
HI and ZEUS should reveal. 

Figure 5 compares our BFKL Born term contributions and the associated re- 
solved contribution with the total NLO cross section da/dxBj- These contributions 
represent more than two thirds of the total NLO corrections in the small XBj region. 
Actually the BFKL Born result of Fig. 5 also contains contributions of graphs in 
which, for instance, the outgoing gluon is attached to the quark line. However these 
contributions are expected to be small as they do not possess the t-singularities (in 
a physical gauge) associated with the exchange of a gluon. To check this point, we 
calculated the contributions of the graph shown in Fig. 4b. For a forward trigger on 
the quark q, we obtain a cross section seven times larger than the one corresponding 
to a trigger on the quark q' (for xsj ~ 2.10 -4 ). Therefore we estimate that the part 
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Figure 6: The BFKL Born term contribution and the total NLO cross section with 
the scales /z 2 = M 2 = Mp = pj_ 4 . Also shown the total NLO cross section with 
scales (Q 2 + E 2 ±4 )/2 (dotted line) and the HI data (4.5 < Q 2 < 15.0 GeV 2 ). 

of the curve of Fig. 5 corresponding to the BFKL Born term of Fig. 4 is dominant. 

From these results we conclude that the main part of the forward cross section is 
due to the BFKL Born terms. Although obtained in the course of a NLO calculation, 
these terms of order 0(a 2 ) represent the Born terms of new channels, namely the 
7*g — > q'q'q and 7*g — > qqg channels. As for any Born terms, we do expect the 
contributions of those channels to be strongly dependent on the renormalization 
and factorization scales. Therefore, contrary to our expectations, we are not able 
to obtain, through our NLO calculations, a total cross section displaying a weak 
dependence on the renormalization and factorization scales. 

Let us be more explicit by studying the effect of the BFKL resummation of the 
small log — terms which appears when extra gluons are emitted by the t-channel 



- BFKL Born terms 

- Total NLO 
Total NLO with scales 
(Q 2 + (dotted line) 
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gluon of Fig. 4 (the BFKL ladder). These contributions have been estimated in ref. 
[SI Ej; they lead to an enhancement of the BFKL Born cross section by a factor 
5 to 10, in obvious contradiction with the data. However the calculations of ref. 
[HI Ej depend on various cuts and do not include the effect of HO corrections to 
the leading BFKL results. These corrections are known to be large; this makes the 
leading results not reliable. A more recent approach |Hj includes a part of these 
HO corrections in its predictions and finds an agreement with HI data [2]. This 
last result, obtained with the scales \l — M — Mp = p±^ (p^ is the momentum of 
the parton which fragments into the 7r°), allows us to make a connection with their 
approach. Using the same scales we obtain the result of Fig. 6 for the BFKL Born 
term and for the total cross section. It is obvious that there is room for a BFKL- 
ladder contribution [H] between the data and the present theoretical prediction. 

However the scale p_i_ 4 is quite large (< p±4 >~ 11 GeV in the HI kinemat- 
ics) compared to what is usually used in large-i?j_ reactions. For instance, in 
the case of 7r° hadropro duct ion in fixed target experiments, a scale M ~ E±±/2 
(< Ej_4 >~ 3.6 GeV in the HI kinematics) is used in ref. [21] to get an agreement 
between data and theory. If a similar scale were used here, we would obtain better 
agreement between data and NLO calculations without any other contributions, as 
is demonstrated in Fig. 6 for the choice fi 2 = M 2 = Mp = \(Q 2 + E\ A ). It is 
clear from this discussion that we cannot accurately determine the importance of 
a BFKL-ladder component in HI data without calculating NLO correction to the 
BFKL Born terms considered here or, in other words, without calculating NNLO 
correction to the electroproduction cross section. 

Conclusions similar to the ones of this section have been obtained by Kramer 
and Potter in their study of the forward leptoproduction of jets l'2oj . 
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Figure 7: The laige-E^ 7r° cross section integrated over E± with the cuts E± > 
3 GeV/c (left) and E ± > 7 GeV/c (right). 

5 Cross section at central rapidity 

In this section we study the leptopro duct ion of large- E±_ ir° in the central region 
in rapidity (in the laboratory frame) with the aim to reduce the BFKL Born term 
contributions and, consequently, to have a better control of the HO corrections. 
Therefore the experimental results obtained in this kinematical domain should pro- 
vide good tests of QCD and the possibility of measuring the virtual photon structure 
function. 

Let us consider the following kinematical range which has been explored by the 
HI collaboration in its measurement of the photoproduction of large- E± hadrons 
namely \JS ep = 300 GeV, .3 < y < .7 and —1 < 7] hadron < 1. We study two 
lower limits for the hadron transverse energy, 3 GeV < E± and 7 GeV < E±, in 
order, as in ref. |2Zj, to estimate the importance of these cuts on the control of the 
HO corrections. For Q 2 , we choose 5 GeV 2 < Q 2 < 10 GeV 2 which belongs to the 
range studied by HI in the measurement of the dijet cross section at low Q 2 fZE\ . 

In Fig. 7 we display the cross section da / di] hadron integrated over E± with E± > 
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Figure 8: The cross section da/dE± for three different choices of scales CJQ 2 + E 



3 GeV (left) and E± > 7 GeV (right). All the scales are set equal to ^Q 2 + E\. 
We clearly see the decrease of the ratio r = -If^ 2 - when the cut on E± increases, 

J Horn - 1 - ' 

with a value r ~ .65 obtained for E± > 7 GeV. On the other hand the ratio of 
the resolved contribution to the Born contributions increases because of the larger 
value of (Q 2 + Ej_)/Q 2 and reaches a value close to 0.5. One must notice that this 
ratio is much smaller than in photoproduction [37] , since the virtuality of the photon 
suppresses the resolved contribution. 

We study the scale variation of the cross-section da/dE± in Fig. 8. Three pre- 



dictions, obtained with the scales ft = M = M F = CyQ 2 + E\ and C — |, 1, 2, are 
displayed. A change of the scales by a factor 4 results in a change of the cross section 
by a factor 2. Therefore, even at large value of E±, the cross section is very sensitive 
to the scale variation. A similar behavior has been observed in the photoproduction 
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cross section of large- E± hadrons [23 • 

6 Conclusion 

In this work we have carried out the calculation of HO corrections to the lep- 
toproduction cross section of large- E± hadrons. These corrections are implemented 
in a parton event-generator offering greater flexibility for the estimation of various 
observables. Concerning the numerical importance of these corrections, we have 
focussed on two different kinematical domains, namely the central region in rapidity 
(in the HERA laboratory frame), and the forward region where we compare our 
results with HI data. 

In the central region we have found important HO corrections at low E±. How- 
ever these corrections decrease as E± increases with a X-factor (HO/ (HO + Born)) 
of about .5 obtained for E± > 7 GeV/c. In the same range the HO corrections 
contain a non-negligible resolved contribution which should allow experiment to 
constrain the virtual photon structure function. 

In the forward region, we find very large HO corrections due to the opening 
of new channels related to the BFKL Born terms ; two thirds of the NLO cross 
section is due to these contributions. These Born terms, and consequently the 
total NLO cross section, are quite sensitive to scale variations and this forbids any 
absolute normalization of the cross section. However one must keep in mind that a 
good agreement between the HI data and the NLO cross section is obtained with 
renormalization and factorization scales taken equal to J (Q 2 + E 2 _)/2. 

Because of the scale sensitivity, no firm conclusion can be drawn on the impor- 
tance of the BFKL resummation in the forward cross section as measured by the HI 
collaboration. Clearly such a study requires the calculation of NNLO corrections to 
the leptoproduction cross section. 
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